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Abstract 

Top quark pair production at proton-antiproton colliders is known to exhibit a forward- backward 
asymmetry due to higher-order QCD effects. We explore how this asymmetry might be studied 
at the Fermilab Tevatron, including how the asymmetry depends on the kinematics of extra hard 
partons. We consider results for top quark pair events with one and two additional hard jets. 
We further note that a similar asymmetry, correlated with the presence of jets, arises in specific 
models for parton showers in Monte Carlo simulations. We conclude that the measurement of this 
asymmetry at the Tevatron will be challenging, but important both for our understanding of QCD 
and for our efforts to model it. 
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I. INTRODUCTION 



The matter content of the highly successful Standard Model (SM) of particle physics was 
generally considered to be fully revealed after the discovery of the top quark in 1994 
While the exact mechanism of electroweak symmetry breaking (EWSB) remains undeter- 
mined 0, the unusually large top quark mass 0, 0], of the same size as the EWSB scale, 
suggests that the top quark may play a special role. Taking a closer experimental look at the 
top quark is therefore high on our list of priorities - not to mention that its quantum num- 
bers (other than mass) and couplings to other SM particles are only crudely known, or not 
known at all [||. Furthermore, the production and subsequent decay of the top quark at cur- 
rent and future collider experiments is a serious background to many new physics searches. 
Thus it is essential to study top quark production in much more detail, and eventually to 
be able to simulate it accurately in Monte Carlo background studies. 

The top quark sample that exists today, exclusively from Runs I and II of the Fermilab 
Tevatron, is rather small, a few hundred events. This shall soon increase by an order of 
magnitude or more. The Tevatron experiments CDF and D0 are beginning to study the 
amount of additional jet activity in top quark pair events, and early indications are that it 
is non-trivial p. There are also hints in the data from Run I that the kinematic distribution 
of the top quark pairs' decay products exhibit a charge asymmetry 0]. 

We reexamine QCD top quark pair production in pp collisions at Tevatron Run II energy, 
with special attention to the asymmetries in the production process and to the correlation 
of the asymmetries with extra radiation in top quark pair events. In fact, the presence of 
both the asymmetries and extra radiation have been theoretically expected for some time, 
but neither the structure of the asymmetries nor their correlations with extra radiation are 
fully reproduced by the standard tools used in experimental analysis. We make the case 
here that it is now time to study the asymmetry and its correlation with the extra radiation 
in detail. In Sec. |H]we review inclusive QCD top quark pair production at hadron colliders, 
including known next-to-leading order (NLO) effects. Previous work had identified an overall 
asymmetry for top quark pairs at leading order (LO) 8] and NLO ;9|, although it did not 
identify all the kinematic features present. We present updated results for the asymmetry 
at NLO, for both inclusive and exclusive samples. In Sec. II 111 we present new results for 
the real emission component of inclusive ti production, where the asymmetry is shown as 
a differential distribution as a function of the kinematics of the extra hard radiation. We 
provide an analysis of the likely statistical uncertainties in a variety of scenarios with various 
luminosities, kinematic cuts and tagging efficiencies. We do not consider the questions 
of systematics and background rates in any detail. In Sections IIVI and IS we study the 
asymmetry for ti production in the zero jet exclusive sample (additional jet activity is 
vetoed) and double jet sample, respectively. Section EU provides a discussion of the level 
of asymmetry that can be found in parton-shower Monte Carlo (PSMC) results where the 
NLO perturbative (coherent) asymmetry effects, as described in the previous sections, are 
not present. Finally we review our results and conclude in Section IV11I 



II. INCLUSIVE QCD TOP QUARK PAIR PRODUCTION 

At LO, ti production is totally charge-conjugation symmetric for both production mech- 
anisms, quark- and gluon-fusion: qq —>■ ti and gg — > ti. As a consequence, the angular 
distributions of the t and i are totally symmetric for pp collisions. However, at higher orders 
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in a s , this no longer remains true for all subprocesses. 

As was pointed out almost two decades ago [8!], not all processes involving additional 
partons are symmetric under charge conjugation with respect to the incoming proton and 
anti-proton beams. The process gg — > tig is, but the processes qq — > tig and qg — > tiq 
are not. Processes involving initial state valence quarks will therefore exhibit a forward- 
backward asymmetry. This is caused by interference between initial- and final-state gluon 
emission (or its crossed process), analogous to what happens in QED, e .g., the forward- 
backward asymmetry in e + e~ — ► /i + yU~7 or other heavy fermion pairs [l0|. Because it 
production at the Tevatron is dominated at the 90% level by qq annihilation, we can expect 
the qq subprocess' asymmetry to be closely reflected in the total sample. 

When the pp — > it cross section was calculated more fully at NLO ^H, an asymmetry 
due to virtual corrections was noticed. It arises from an interference between the color- 
singlet 4-point (box) virtual correction and the Born term for the process qq — > ti. Ref. [j| 
examined this asymmetry more closely, although still inclusively and integrated over the 
phase space of the additional parton. It was found that the virtual contribution produces 
an asymmetry opposite in sign to and larger than that of the real emission component. For 
the forward-backward inclusive asymmetry defined in terms of the top quark rapidity yt l , 

N t (y t >0)-N t (yt < 0) 
FB N t (y t >0) + N t (y t <0) ' {) 

Ref. ^ calculated a value of 4 — 5% for the Tevatron Run I, y/s = 1.8 TeV, and top quark 
mass m t = 175 GeV. 



A. Recalculation of top quark asymmetry 

To provide a connection to the previous work outlined above, we have recalculated the top 
and anti-top quark distributions for Tevatron Run II, \/s = 1.96 TeV, using the NLO Monte 
Carlo code mcfm with the updated top quark mass, m t = 178 GeV, and structure 
functions of CTEQ6 [13| (set LI for LO, M for NLO). We show our results in Fig. ^for the 
LO and NLO inclusive calculations, and also for the tiOj exclusive calculations. We define 
this last quantity as the NLO ti inclusive rate minus the LO tij inclusive rate above some 
cutoff, ptU) > Prijifnin) and inside some rapidity region \r](j)\ < r](j,max). (Note that 
inclusive NLO ti and LO tij are the same order in a s .) Thus we have 

°m] = 4- L ° - 4°(ptU) > PT(j,min), \ V (j)\ < r,(j,max)) . (2) 

The tiOj exclusive rate corresponds to the total ti + X rate with a jet veto applied above 
some pt inside some fiducial region of rapidity in the detector. Here we use the cutoffs 

Mj)>20GeV, \ V (j)\<3.0. 

This fiducial region is based on the Run II capabilities of the CDF and D0 detectors. It 
is important to note that, in order to evaluate the perturbative tij asymmetry described 



Here we use rapidity, y, instead of polar angle, 9, since rapidity is conventional for hadron collider exper- 
iments. For massless particles, rapidity y is equal to pseudorapidity i]. 
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above, the calculation must employ the exact matrix elements. Since parton-shower Monte 
Carlo tools such as pythia lj, herwig l]| and sherpa [l6| do not include the full 



interference effects in the parton radiation that produce the asymmetry, they cannot be 
used to predict the asymmetry. However, as we will discuss later, the soft radiation models 
used in the Monte Carlos are capable of producing certain kinds of asymmetries, which 
must be understood in order to accurately analyze any asymmetries that may be in the 
experimental data. In our perturbative calculations the minimum p? kinematic cut protects 
us from the soft singularity present in the tig final state. We verify that this is sufficient 
by checking that a t tj <C <Jtt- For the kinematic cuts defined above, mcfm gives a t i = 5.9 pb 
and a t tj = 1.1 pb. Their ratio is only 0.19, indicating that the calculation is perturbative. 
However, because a t tj is calculated only at LO, there is a considerable residual uncertainty of 
probably a factor of two at the Tevatron energy. This means there is also a larger uncertainty 
in the tiOj exclusive rate compared to the NLO inclusive rate, since the former is effectively 
calculated at LO. For example, we choose the top quark mass, m t , as the factorization and 
renormalization scales in our benchmark cross section. This is conservative because it lies 
at the low end of our estimates for the LO tij rate. However, if we choose scales that take 
into account the softness of the additional radiation parton, e.g., pr(j), the cross section 
can be twice as large as our benchmark point. Even in this case we would still regard our 
result as perturbative. For a final analysis from Run II data, knowing the NLO tij rate will 
probably become necessary. Fortunately, this calculation is underway ET 
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FIG. 1: Differential cross section distributions as a function of the top (left) and anti-top (right) 
quark rapidities, produced in pp collisions at the Tevatron in Run II, \fs = 1.96 TeV. Shown are the 
LO ti inclusive (dotted), ti NLO inclusive (solid), LO tij inclusive (dot-dashed) and tiOj exclusive 
(dashed) predictions. We define the LO tij inclusive rate as that where the additional final-state 
parton has ptU) > 20 GeV and \rj(j)\ < 3. The tiOj exclusive rate is then the NLO inclusive rate 
minus the LO tij inclusive rate. 
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The structure of the asymmetry described above can be seen in Fig. ^ First note that 
the LO ti inclusive cross section (dotted curve) is symmetric about y = in both panels, i.e., 
for both t and i. Comparison of this symmetric curve with the (solid curve) NLO inclusive 
result indicates that the latter curve is (slightly) shifted to larger y for t (left panel) and 
smaller y for i (right panel) corresponding to a positive asymmetry. The tiOj exclusive cross 
section (dashed curve) exhibits a shift that is similar in terms of magnitude and direction. 
On the other hand, the LO tij inclusive cross section (dot-dashed curve) is shifted in the 
opposite direction in each panel yielding a negative asymmetry. Thus we conclude that the 
real emission corrections tend to push the top quark backward, opposite to the proton beam 
direction, while the virtual corrections push it forward, in the direction of the proton beam. 
The virtual corrections are larger as indicated by the qualitative agreement between the 
overall NLO inclusive distribution (solid curve) and the tiOj exclusive distribution (dashed 
curve). The y t and yi distributions are mirror images of each other, as required by CP- 
invariance of QCD. 



process 


A 1 
FB 


LO ti inclusive 





NLO ti inclusive 


3.8 


LO tij inclusive 


-6.9 


tiOj exclusive 


6.4 



TABLE I: Top quark forward-backward asymmetry in pp — > ti production at Tevatron Run II, 
= 1.96 TeV, for top quark mass m t = 178 GeV. The lj inclusive and tiOj exclusive rates are 
defined by pr(j) > 20GeV and \r](j)\ < 3.0. A FB = —A FB by CP-invariance. The cross sections 
were calculated with mcfm 

Table U shows our numerical results for A FB that arise from integrating over the dis- 
tributions in Fig. Note that, as expected, these estimates of the magnitude of the tt 
asymmetries depend on the choice of factorization and renormalization scales. For example, 
with both scales set to m t /2, mt, and 2mt, the NLO inclusive ti asymmetries are 5.5%, 
3.8%, 2.8%, respectively. In this paper, we examine all asymmetries with both scales set to 
m t . Since the strong interaction is separately C and P symmetric and the initial state in 
the pp process is a CP eigenstate, overall CP -invariance requires that A FB — — A FB . At 
the same time, since the initial state is not an eigenstate of C (or P) alone, we can have 
A FB ^ A FB 7^ 0. As expected from the distributions, the LO asymmetry vanishes, the 
NLO inclusive and tiOj exclusive asymmetries are positive, and the asymmetry result for 
the tij inclusive distribution is negative. When considering these results, it is important 
to recall that the magnitudes (but not the signs) of the last two entries in the table depend 
on the pt cut defining the jets, which in this case has the value pr{j) > 20 GeV. NLO final 
states with extra real emission above the cut are included in the LO tij inclusive sample 
and exhibit a negative asymmetry. The corresponding real emission configurations below 
the cut are included in the tiOj exclusive result, where they partially cancel the virtual 
contributions with a positive asymmetry reducing the net asymmetry in this term. Thus 
we can vary the magnitudes of the last two entries in the table up and down together by 
varying the jet-defining px cut (larger magnitudes for a smaller cut and vice versa). Our 
numerical results differ slightly from those of Ref. due partly to the higher energy, but 
also due to the different PDF set used here. 
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The obvious question to ask is, can Tevatron Run II measure the various inclusive and 
exclusive asymmetries predicted by perturbative QCD? Observing the top quark asymmetry 
in a collider experiment such as Tevatron Run II is complicated by the difficulty in accurately 
reconstructing the underlying parton kinematics from the measurements of jets, compared 
to leptons. Because top quarks are produced at the Tevatron with transverse momentum 
typically smaller than their mass, their decay products are distributed over a large solid 
angle. Mismeasurement of individual jets thus results in a non-negligible uncertainty in the 
top quark direction. One should then look instead at the decay products, since any forward- 
backward top quark asymmetry will manifest itself in the daughter particles as well, albeit 
with possibly different magnitudes. The b jets are tagged, but their charge is extremely 
difficult to determine, so they are not good candidates. Nor are hadronically-decaying W 
bosons, again because of jet mismeasurement issues, but also because light-flavor jet charges 
cannot be determined at all. Instead we propose to use the asymmetry of the final state 
leptons (electron or muon) from W decays, since the the direction of leptons is extremely 
well-measured. 

Top quark pairs decay to two possible final states containing a lepton: both top quarks 
can decay to b£u, where i = e, /i, which we call the "dilepton" final state; or one decays to 
biv while the other decays to bjj, which we call the "lepton +jets" final state. The cleanest 
sample experimentally is the dilepton sample, although it occurs at only about 1/6 the 
rate of the lepton+jets channel. Assuming both b jets are tagged, the dilepton sample has 
practically no background from other SM processes, while the lepton+jets channel has a 
signal to background (S:B) ratio of about 6:1 [l8j]. The dilepton channel has the added 
advantage that there is no confusion between the additional radiated hard jet in the event 
and jets from the W boson decay, as in the lepton+jets sample. Since these events include 
jets along with the leptons, we will find it informative to examine the correlations between 
the lepton forward-backward asymmetry and the kinematic properties of the jets. 

To obtain a handle on the expected lepton asymmetries we calculate using matrix elements 
for pp — ► ti production, including decay spin correlations and treating the intermediate W 
bosons off-shell, as generated by madgraph [lj|. We evaluate cross sections at \fs = 
2.0 TeV and use factorization and renormalization scale choices of /i/ = fi r = m t (the small 
dependence on the total energy and the factorization/renormalization scales will not be an 
issue here). Before top quark decays, the cross section is 4.87 pb, compared to 6.7-8.0 pb 
for the NLO+(N)NLL ti rate [20]. We will normalize our results to the NLO rate for our 
inclusive results; using the average of the two NLO results, 7.4 pb, the effective K-factor is 
1.52. 

To calculate cross sections for observable final states at the parton level, we use our LO 
madgraph code normalized to the NLO rate, then impose generic kinematic cuts on the 
outgoing particles suitable for both DO and CDF: 

j) T > 20 GeV , 

p T (6)>15GeV, |77(6)|<2.0, AR(b, b) > 0.4 , 

p T (j)>20GeV, \ V (j)\ < 2.0, AR(j,j;j,b)>0A, (3) 

p T (£)> 20 GeV, \r](£)\ < 1.1 , AR(£, £) > 0.2 , AR(£,j;£, b) > 0.4 . 

The cuts are very conservative compared to the CDF and DO expectations for the majority 
of Run II, but reflect the present level of understanding of the detectors and are used in 
current top quark analyses 0, 2^]. For example, muons can easily be identified out to 



a pseudorapidity of 2.0, and non-6 jets can be identified for pr(j) > 15 GeV. For both 
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the dilepton and lepton+jets final states we discuss below, we consider the possibility of 
increasing the statistics by loosening the cuts accordingly. We do not attempt to simulate 
detector effects. 

Regardless of the decay mode, we define the lepton asymmetry as 

t N( Ve >0)-N( Ve <0) 

FB N(r ]e >0) + N(r ]e <0)> {) 

where the lepton may be £ + or £~, depending on which is visible (or both, in dilepton events, 
but A l p B and A FB are calculated separately). Once a cross section with cuts is obtained, 
the NLO-normalized number of expected events for a given luminosity and set of detector 
ID efficiencies is distributed into forward and backward bins such that the asymmetry is the 
same 3.8% as that predicted at NLO by mcfm for the top quarks. It is important to point 
out, as we will see later explicitly in the tij Section ITTT] that kinematic cuts can alter the 
asymmetry for a given subsample. Hence our procedure is only an approximation. However, 
as no NLO tt program yet has the capability to decay the top quarks with spin correlations 
intact at NLO, our approximation is the best we can accomplish at the moment. Hopefully 
it will be superseded in the near future. 



B. Dilepton final state 

When both top quarks decay leptonically, the final state in the detector is bb£ + £~^>T, 
where £ = e, fi and fa is missing energy transverse to the beam axis, resulting from the 
neutrinos which escape undetected. We assume that both b jets are tagged, rendering the 
backgrounds negligible. This is thus the cleanest top quark pair sample in the data. The 
cross section with dilepton kinematic cuts and NLO K-f actor applied is 113 fb. 

We begin with a very conservative baseline assumption of 20% for the double 6-tagging 
efficiency, which we label 62b- For 4 fb -1 and straightforwardly combining the statistics of 
D0 and CDF (i.e., a factor of 2), we estimate about 180 events in the dilepton sample with 
two b tags summed over £ = e, \i . The statistical uncertainty on an asymmetry is given in 
Ref. @: 




where N F (N B ) is the number of observed events with the lepton £ forward (backward). 
Taking CP as a good symmetry we have A l FB = —A FB (N F + = N B , N l B = N F ) and 
we can combine the statistics of the £ + and £~ channels. For the dilepton sample this 
means that each event contributes twice. Combining the two charges also yields a factor 
\[2 improvement in overall significance of the asymmetry determination over the single sign 
analysis. We summarize the dilepton channel results in Table |H] We see that the initial 
numbers for the dilepton channel yield an uncertainty in the asymmetry measurement of 
6A e FB /A e FB ~ 1.4, providing only a 0.7a effect with the summed sign data. We conclude 
that it will be challenging to use the dilepton channel due to poor statistics, at least with our 
conservative assumptions. However, there are many opportunities for increasing the sample 
size. 

Here we explore various prospects for improving the sample size coming from multiple 
sources, of varying likelihood. The results are illustrated in Table |TT| First, the lepton 
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cuts can be loosened by observing leptons out to a pseudorapidity of 2.0 instead of 1.1, 
which is already possible for muons. This adds ~ 65% to the observed rate. Lowering the 
Pt cut of only one lepton to 10 GeV adds another 15%. Another few percent improvement 
will automatically come from events where one or both W bosons decay to a tau lepton, 
which then decays to an electron or muon. Imposing looser lepton cuts but ignoring the 
contribution from taus, we obtain a rate enhancement of about a factor of 1.9 so that the 
uncertainty decreases slightly, to 5A FB /A FB ~ 1.0. The measurement becomes 1 cr as 
indicated in the second line of Table |HJ 

We also expect that the 6-tag efficiency should be improved noticeably, both by adding 
kinematic track information in a neural net analysis to the vertex tag, and by beginning 



to include soft lepton tags, which is about 10% per b [2l|. Together, these highly likely 
improvements to fe-jet tagging would enlarge the sample by another factor of 1.9. Third, 
Run II could achieve its "stretch" goal of about 8 fb _1 , which would double the sample. 
Taken together (the most optimistic scenario), these improvements would result in a dilepton 
sample more about seven times what our conservative scenario predicts. In this case, the 
asymmetry in the dilepton channel is still challenging, yielding a 2.0 a measurement of the 
inclusive asymmetry as indicated in the fourth line of Table |H] 

It is fair to ask how feasible these improvements are. The expanded lepton pseudorapidity 
coverage is a safe bet, while the lowered- pr cut is possible but not guaranteed (and also 
doesn't result in much gain, compared to increased angular coverage). The 6-jet tagging 
efficiency will definitely improve noticeably with the addition of soft lepton tags, while the 
prospects for improved silicon vertex tags are not yet fully understood. We cannot guess at 
the total integrated luminosity of Run II, but anticipate continued efforts at improvement. 

Another natural possibility to consider is requiring only one 6-tag, which with the assumed 
improved tagging efficiency would increase the sample by about a factor of 2.2. While this 
approach introduces some SM background, the background contribution will not be large. 
Ultimately this must be taken into account, although it is beyond the scope of this paper. 
Ignoring the small contribution from background, a single 6-tag strategy would result in 
an uncertainty of 5A FB /A FB ~ 0.34 and yield a 2.9 a measurement as indicated in the 
fifth line in Table El We note that it may also be possible to usefully study a zero 6-tag 
dilepton sample. Half this sample, containing an e + /i mixed-flavor pair, has extremely low 
backgrounds. The primary background for the other half, ee or /z/i events, arises from the 
Drell-Yan process with an intermediate Z boson with a substantial (~ 20%) asymmetry. We 
expect this background can be controlled by cutting away the Z-pole region and demanding 
extra jets. 

We also mention one other possible scheme for improving the measurement in the dilepton 
sample. Instead of an asymmetry based on a single lepton being forward versus background 
in the laboratory frame, consider one based on the longitudinal ordering of the 2 leptons, 

A i N(yi+ > Vi-) - N fa+ < rie-) 
FB N( Vi+ > m -) + N(r ]e+ < Ve -) ' 1 ' 

A preliminary look at the inclusive asymmetry using mcfm [l^] suggests that A l FB may be 
50% larger than A l FB . While the dilepton sample is likely not as statistically significant as 
lepton+jets, this observable may become more relevant if systematic uncertainties prove to 
be the larger problem. 
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Sample 


JCdt [fb- 1 ] 


cuts 


£26 


< + n£ 


Ng + Nf 


SA FB 
A FR 


signif. 


dilepton 
dilepton 
dilepton 
dilepton 
dilepton 


4 
4 
4 
8 
8 


tight 
loose 
loose 
loose 
loose 


0.20 
0.20 
0.38 
0.38 
0.62* 


188 
360 
685 
1370 
3074 


174 

334 

635 

1270 

2849 


1.4 
1.0 
0.72 
0.51 
0.34 


0.7 a 
1.0 a 
1.4 a 
2.0 a 
2.9 a 


lepton+jets 
lepton+jets 
lepton+jets 
lepton+jets 


4 
4 
4 

8 


tight 
loose 
loose 
loose 


0.20 
0.20 
0.38 
0.38 


570 
864 
1642 
3285 


528 
801 
1522 
3044 


0.79 
0.64 
0.47 
0.33 


1.3 a 
1.6 a 
2.1 a 
3.0 a 


combined 
combined 
combined 
combined 
combined 


4 
4 
4 
8 
8 


tight 
loose 
loose 
loose 
loose 


0.20 
0.20 
0.38 
0.38 
0.38* 


757 
1225 
2327 
4654 
6359 


702 

1135 

2157 

4314 

5893 


0.69 
0.54 
0.39 
0.28 
0.24 


1.5 a 
1.8 a 

2.5 a 

3.6 a 
4.2 a 



TABLE II: Numbers of forward and backward lepton events, combining £ + and l~ samples (rounded 
to the nearest integer) and expected statistical uncertainty on the measured lepton asymmetry, 
A FB , for inclusive ti production in pp collisions at Tevatron Run II, y/s = 2.0 TeV, for top quark 
mass mt = 178 GeV, summed over the two detectors. The lepton asymmetry in each case is 
assumed to be that of the top quarks themselves, +3.8%. The upper block contains the result 
for the dilepton sample, the middle block for the lepton+jets sample, and the lower block for the 
combined samples. The total NLO cross section is taken to be 7.4 pb, the average of the results of 
Refs. [20j. "Tight" cuts refers to those of Eq. [21 while "loose" refers to the increased acceptance 
scenario described in Sees. Ill Bill Cl The * entries represent using a single 6-tag strategy for the 
dilepton sample. 

C. Lepton+jets final state 

A different approach to increasing the sample size is to consider the lepton+jets channel, 
which is slightly more than a factor six larger in branching ratio than the dilepton channel. 
After cuts, the LO cross section times NLO K-factor is 686 fb, approximately maintaining 
the ratio of lepton+jets to dilepton rates. This channel does have a larger background 
than dileptons, mostly from VT+jets, which is not fully understood either theoretically or 
from data. However, with double 6-jet tagging, S:B is about 6:1 [l8j], so for our purposes 
we assume a double 6-tagging strategy and ignore background contamination. Naturally, 
the inherent asymmetries of the same final state from VT+jets, about 20% in the untagged 
sample, must eventually be included, using exact matrix elements also for the background. 

Again conservatively assuming a double 6-jet tagging efficiency of 20%, we predict about 
1100 events in 4 fb -1 summed over both experiments. Note that for this channel each event 
only contributes once (a single lepton) to the measurement of the asymmetry, for a net 
improvement over the dilepton case of about a factor of 3. As indicated in the middle 
block of Table HP this provides for an initial measurement: 5A FB /A FB ~ 0.8 and statistical 
significance of 1.3 a. As above we consider the impact of increasing the lepton acceptance 
range to r)(£) < 2.0, lowering the lepton pr(J) cut to 15 GeV, and allowing for non-6 jets 
from W decay also down to pt > 15 GeV. This would increase the cross section by about 
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50% as indicated in the middle block of Table ITT1 Prospects for improving this sample's size 
from other sources, such as improved 6-tagging and greater integrated luminosity, are the 
same as in the dilepton sample. In the most optimistic scenario of 8 fb _1 per experiment, a 
double 6-tagging efficiency of 38% and utilizing the looser cuts, we find that the uncertainty 
in the asymmetry would be about 1/3, leading to a 3a measurement as indicated in Table ITT1 
Combining this and the dilepton channels would increase the statistical significance to more 
than 4<j as indicated in the bottom section of Table |H] 



III. DIFFERENTIAL ASYMMETRY IN TOP QUARK PAIR +1 JET EVENTS 

Next we consider correlations between the lepton asymmetry and a single extra jet in the 
event not arising from the top quark decays. As in the inclusive case, for all our calculations 
we use matrix elements for pp — > tij production (j = q,q,g summed over four light quark 
flavors), including decay spin correlations and treating the intermediate W bosons off-shell, 
as generated by madgraph The extra parton does not come from either of the on- 
shell top quark decays, i.e., we do not include radiative decays. We calculate cross sections 
at ,/s = 2.0 TeV requiring the additional hard jet to have px > 20 GeV and \r]\ < 2.0, 
and use factorization and renormalization scale choices of fif = fi r = m t , which yields a 
perturbatively well-behaved result: before top quark decays, the cross section is 1.03 pb, 
compared to 6.7-8.0 pb (avg. 7.4 pb) for the NLO+(N)NLL ti rate @. 

As noted in the previous section, choosing a more "physical" scale which takes into 
account the softness of the additional parton, the cross section rises to 2.0 pb, still a rea- 
sonable, perturbative result. The uncertainty on the inclusive tij rate awaits a full NLO 
calculation 17], but our conservative choice of /i = m t gives the lowest cross section and 
may be regarded as a conservative baseline. We will consider the possibility of larger nor- 
malization (by a factor of 2) in our numerical results. Note that our predictions for the 
asymmetry will not depend on the normalization, as discussed in the next section. 



A. Dilepton final state 

We assume that both b jets are tagged, so it is trivial to identify the additional parton. 
The additional parton in the data, however, will not always arise from the tij "production- 
radiation" process as we have defined it. Approximately 50% of the time, i.e., at a rate 
approximately equal to what we calculate here, the extra jet will be due to radiative top 
quark decay |24|. In the dilepton sample, this occurs when one of the b quarks emits a hard 
gluon. This defines the "decay-radiation" contribution. (As argued in Ref. 24]| this separa- 



tion between production- and decay-radiation, while approximate, is numerically reliable.) 
We discuss later how the production and decay contributions might be separated. As the 
calculations here describe the appearance of extra radiation in the production mechanism, 
they will exhibit the asymmetry characteristic of the LO tij inclusive sample. 

We could again show the overall asymmetry, as in Sec. EJ integrated over all kinematics 
of the additional jet, but this is just a single number, with magnitude about 4 — 5%, the 
same as the asymmetry of the top quarks themselves. We do note that when the top quarks 
decay, leptons can be thrown preferentially back into the hemisphere opposite to the top 
quark flight direction, due to the spin of the top quark. Thus, the lepton asymmetry will 
not necessarily be the same as the top quark asymmetry show in Sec. |H] before decays. In 
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any case, the more interesting result lies in the asymmetry as a function of the kinematics 
of the additional jet. We first show A FB as a function of the pr of the extra jet in Fig. |21 
The corresponding version of Eq. 0] 



A FB (pt {])) 



Njrn > 0,p T (j)) -N{ m < 0,p T {j)) 
N{r) £ > Q,p T (])) + N (rj t < 0,p T (])) 



(?) 



where, for example, N(rj£ > 0,pr(j)) is the number of events with a forward lepton and 
an extra jet with transverse momentum px (j). As indicated in Fig. this distribution is 
essentially constant at the overall asymmetry value. Despite the lack of structure, this result 
is important because it shows that the asymmetry is not dependent on the extra parton px 
cutoff, so our results are independent of uncertainty due to the choice of the perturbative 
cutoff in p T . CP still tells us that A FB (p T (j)) = —A FB (p T (j)) as shown in the figure. 
Note also that, as expected, A e FB < 0, just as for the LO tij inclusive sample. 
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FIG. 2: Forward-backward lepton asymmetry in production-radiation tij dilepton events as a 
function of the transverse momentum of the additional hard jet. The £ + (£~) distribution is the 
solid (dashed) curve. The two curves are CP-invariant up to the level of Monte Carlo statistical 
uncertainty. 



The asymmetry with respect to the angular distribution of the additional jet is far more 
revealing. The distribution of interest is now 



N(ru>Q,ri(j))-N(ri e <0,ri(j)) 
N(r le >0,r l (j)) + N(r le <0,r l (j)) 



(8) 



arc 



where CP invariance specifies that A FB (j](j)) = — A FB (— r\ (j)), i.e., A e FB and A l FB 
mirror images of each other about r](j) = 0. We show this asymmetry distribution in 
Fig. El separately for the different parton subprocess and for the combined result. Note 
that for each pair of distributions (for the 2 lepton charges) the CP relation just noted 
is satisfied, but not always in the same way. For the subprocess gg — > tig the initial 
state is symmetric with respect to C alone requiring that A e FB (rj (j) , gg) = A FB (77 (j) , gg), 
so that CP now requires that the distributions are separately odd in the jet rapidity, 
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FIG. 3: Forward-backward lepton asymmetries for each subprocess in production-radiation ttj 
dilepton events as a function of the pseudorapidity of the additional hard jet. The t + (£~) distri- 
butions are shown by the solid (dashed) curves. The left panel shows the dominant qq contribution 
to the asymmetry (the curves with the slightly larger magnitude for both £ + and i~) and the 
asymmetry for both charges for the total rate (the curves with the smaller magnitude). The small 
difference between qq and total arises from the contributions of the other parton channels., whose 
asymmetries are indicated in the right panel. Note the absence of any charge dependence in the 
curves in the right panel. 



A^ B (77 (j) , gg) = —A FB (—r](j),gg). These constraints are in agreement with the cor- 
responding (middle) curves in the right panel of Fig. |H1 The specific form of these dis- 
tributions suggests that for this subprocess the ti pair tends to recoil from the extra jet, 
i.e., the t and i quarks (and the leptons from their decays) tend to be in the same hemi- 
sphere while the recoiling jet is in the opposite hemisphere (a backward jet yields a pos- 
itive asymmetry and conversely). Formally, the subprocesses containing a single quark 
line exhibit an inclusive asymmetry, as found in Ref. [9], although it is numerically ex- 
tremely small. For the discussion here we note that the combinations qg + qg or gq + gq, 
as in the upper and lower curves in the right panel of Fig. are strictly not C sym- 
metric states due to the differences between q and q distributions within a proton (or 
antiproton). Thus it does not follow that A FB (rj (j) , qg + qg) = A FB (77 (j) ,qg + qg), 

A E fb (v ti) ' 9Q + 5"?) = ^fb (v ti) ) 91 + 90)i even though these relations are approximately 
numerically true in our results in Fig. El So in this case the (exact) constraint of CP leads 
only to the approximate relation A e FB (77 (j) , qg + qg) ~ — A FB {—rj (j) , gq + gq). These fea- 
tures are illustrated in Fig. |21 where the boost effects from the vastly different average Feyn- 
man x values for quarks and gluons in the proton (antiproton) result in the specific shapes of 
the curves. Qualitatively the qg + qg subprocess tends to a produce a ti pair in the forward 
hemisphere (A FB > 0) with the extra jet recoiling into the backward hemisphere. For- 
ward/backward are reversed for the gq + gq subprocess. Note that, after summing the gg, qg 
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and gq subprocesses, their net asymmetry integrates to approximately zero (over a symmetric 
interval in T)(j)). Finally consider the qq subprocess, which not only dominates the total rate 
at Tevatron energies, but also exhibits the most interesting structure. The qq initial state 
is not separately invariant under C and P so that we have A^ B (77 (j) , qq) 7^ A FB (77 (j) , qq), 
while CP still guarantees that A FB (77 (j) , qq) = —A FB (—i](j),qq). The corresponding 
asymmetry distributions are illustrated in the left panel of Fig. El for both the qq subprocess 
(the curves with large magnitudes) and the total of all subprocesses (the curves with slightly 
smaller magnitudes). Note that for these distributions the asymmetry for a single lepton 
charge clearly does not integrate to zero. Qualitatively the qq subprocess tends to produce 
a t quark in the backward hemisphere and a i antiquark at rest in the lab when the extra 
jet is forward, while a backward jet tends to correlate with a t antiquark in the forward 
hemisphere and a t quark at rest in the lab. 

As indicated in Fig. [HI the asymmetry reaches a maximum absolute value of around 14% 
for extra jet rapidity in the region around \t]{j)\ ~ 0.4, and, even averaged over the region 
< r]{j) < 1 (for £ + ), the asymmetry is over 10%. The question of whether this can actually 
be measured is more difficult to answer. The dilepton production-radiation ttj cross section 
at LO with the "tight" cuts of Eq. is 15.1 fb. With 4 fb" 1 and a conservative double 6-tag 
rate of 20%, Run II will then observe about 24 dilepton production-radiation ttj double-6- 
tagged events for D0 and CDF combined, based on the LO rate (and an approximately equal 
number of radiative ttj events). This is not enough statistics to perform even an overall 
asymmetry measurement, much less in a restricted region of i](j). However, with the 
prospects for improving the sample size as discussed in Sec. Ill Bl as well as the uncertainty 
in cross section normalization, this channel could possibly collect (9(100) events and would 
then become interesting. 

We present Fig. the normalized angular-differential cross section distribution overlaid 
on the (positive) lepton asymmetry, as a guide to how binning in 77 (j) might be done depend- 
ing on the ultimate statistics achieved in the dilepton channel. Fairly central jet rapidity 
features the largest asymmetry as well as large portion of the total rate. At a minimum, one 
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FIG. 4: Normalized differential cross section with respect to the extra jet pseudorapidity (dashed), 
overlaid with the total forward-backward positive lepton asymmetry (solid), in dilepton production- 
radiation ttj events with loose cuts as described in the text. 
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could use two bins, jet-forward and jet-backward, discarding the region beyond \rjj\ ~ 1. 

We have avoided the issue of the approximately equal-sized sample of radiative decay 
tij dilepton events. This is obviously a non-trivial complication. Since the statistics of 
the dilepton channel are weak to begin with, we don't attempt to address it further here, 
except to note that radiative decay events are highly suppressed by imposing an angular 
cut ARj b 3> 0.4. For exam ple, the cut ARj b > 1.0 reduces the radiative decay cross section 
by more than a factor two [24], but our calculations show only a 20% loss of production tij 
dilepton events. In the lucky circumstance that Run II collects a very large dilepton top 
quark sample, this issue should be thoroughly investigated. 



B. Lepton+jets final state 

We again also consider the lepton+jets channel, which is slightly more than a factor six 
larger in branching ratio than the dilepton channel, although it does suffer from a slight 
background as discussed in Sec. Ill CI After cuts, the LO cross section is 88.9 fb, slightly less 
than a factor six larger than the dilepton channel with cuts. 

The lepton+jets channel has one complication, that of correctly identifying which of the 
three non-6-tagged jets is the "extra" one, i.e., the one that is not a top quark decay product. 
This situation is aggravated by the combinatorics, but can be largely addressed by imposing 
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FIG. 5: Forward-backward lepton asymmetries for each subprocess in production-radiation tij 
lepton+jets events as a function of the pseudorapidity of the additional hard jet. The £ + (£~) 
distributions are shown by the solid (dashed) curves. The left panel shows the dominant qq 
contribution to the asymmetry (the curves with the slightly larger magnitude for both £ + and t~) 
and the asymmetry for both charges for the total rate (the curves with the smaller magnitude). The 
small difference between qq and total arises from the contributions of the other parton channels., 
whose asymmetries are indicated in the right panel. Note the absence of any charge dependence 
in the curves in the right panel. 
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the W mass constraint on one jet pair, the top quark mass constraint on that pair and a b 
jet, and the transverse top mass constraint on the other b jet, the lepton and the missing 
transverse energy. We do not examine this issue further here, but expect misidentification 
effects due to kinematical combinatorics to be a minor correction. This issue should of course 
be examined in detail when the channel is studied with full detector simulation. Such mass 
constraints also remove the vast bulk of top quark decay-radiation events from the sample, 
as our minimum Pt{j) requirement means that sample contamination from radiative events 
enters mostly via jet mismeasurement at the edges of invariant mass windows. 

Fig.EJshows the lepton asymmetry in production-radiation tij lepton+jets events, similar 
to Fig. El for the dilepton case, also with the cuts of Eq. El This plot assumes that the 
additional jet, which did not come from a top decay, is determined with 100% accuracy, as 
discussed above. The asymmetry features in this sample are almost identical to those in 
the dilepton sample. The asymmetry as a function of extra jet is flat as in the dilepton 
channel, and we do not show the corresponding figure here. 

Again starting from the assumption of 20% efficiency for double 6-jet tagging (no im- 
provement over the current level), our LO estimate predicts about 142 events in 4 fb _1 for 
both experiments. Measuring the differential asymmetry therefore appears possible only for 
large pseudorapidity bins. Applying the loose cuts as in the inclusive sample gives about a 
factor 1.5 increase. Note that we keep the requirement that the extra jet has pt > 20 GeV, 
simply to keep our calculation confidently in the perturbative regime. Prospects for im- 
proving this sample's size from other sources, such as improved 6-tagging, greater integrated 
luminosity, or cross section normalization, are the same as in the dilepton sample. 

C. Numerical results 

Here we make some simple estimates of the uncertainty that might be achieved on an 
A l FB measurement in the inclusive tij sample. The ideal objective in this case would be to 
map out the asymmetry with respect to the angular structure of the additional hard parton, 




FIG. 6: Normalized differential ttj cross section with respect to the extra jet pseudorapidity 
(dashed), overlaid with the total forward-backward positive lepton asymmetry (solid), in lep- 
ton+jets events with loose cuts as described in the text. 
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which is additional information that the inclusive and exclusive tt samples do not have. 
We consider both the dilepton and lepton+jets samples, although the latter has far better 
statistics and is expected to largely avoid the complications of radiative top quark decay. 
We consider only the statistical uncertainty, as the most important systematic uncertainties 
are related to detector effects, which we cannot include; or the very small backgrounds 
in double 6-tagged events, which are expected to be a much smaller uncertainty than the 
limited statistics will allow for. 

Due to the limited statistics available, we select a single bin in extra jet pseudorapidity of 
—0.1 < rj(j) < 1.1 (—1.1 < < 0.1) for £ + (£~). Varying this over a range of a few tenths 
either direction at either limit does very little to the overall uncertainty, while increasing 
the range to increase statistics results in an overall lower asymmetry, and vice versa, as seen 
in Fig. The asymmetry in this bin is 11%. 

We show our numerical results in Table ITTT1 which has a general format similar to Table ITT1 
As suggested by Figs. |U and El approximately 1/3 of the events fall into the rapidity window 
noted above, —0.1 < T)(j) < 1.1 (—1.1 < 77 (j) < 0.1) for £ + (£~). Thus, since the dilepton 
events still contribute twice, the 24 dilepton events mentioned above yield 16 counted leptons 
in the first line of the dilepton block (24 x 1/3 x 2) in Table llTH For the lepton+jets sample the 
number of leptons is just 1/3 of the number of events so that 142 events become 47 leptons. 



Sample 


JCdt [fb- 1 ] 


cuts 


£26 




Np + N£ 


<+ Nf 


K ttj 

NLO 


A WR 


signif. 


dilepton 
dilepton 
dilepton 
dilepton 
dilepton 


4 
4 
4 
8 
8 


tight 
loose 
loose 
loose 
loose 


0.20 
0.20 
0.38 
0.38 
0.38 


-0.11 

-0.12 
-0.12 
-0.12 
-0.12 


7 

14 
26 
51 
103 


9 

17 
33 
66 
132 


1 
1 
1 
1 

2 


2.3 
1.5 
1.1 
0.75 
0.53 


0.4 a 
0.7 a 
0.9 a 
1.3 a 
1.9 a 


lepton+jets 
lepton+jets 
lepton+jets 
lepton+jets 
lepton+jets 


4 
4 
4 
8 
8 


tight 
loose 
loose 
loose 
loose 


0.20 
0.20 
0.38 
0.38 
0.38 


-0.11 
-0.12 
-0.12 
-0.12 
-0.12 


21 

31 

60 

119 

239 


26 

40 

76 

153 

306 


1 
1 
1 
1 

2 


0.8 
1.0 
0.69 
0.49 
0.35 


0.8 (J 
1.0 a 
1.4 a 
2.0 a 
2.9 a 


combined 
combined 
combined 
combined 
combined 


4 
4 
4 
8 
8 


tight 
loose 
loose 
loose 
loose 


0.20 
0.20 
0.38 
0.38 
0.38 


-0.11 
-0.12 
-0.12 
-0.12 
-0.12 


28 
45 
85 
171 
342 


36 

58 

109 

219 

437 


1 
1 
1 
1 

2 


1.1 

0.8 
0.58 
0.41 
0.29 


0.9 a 

1.3 a 
1.7 a 

2.4 a 
3.4 a 



TABLE III: Numbers of forward and backward lepton events combining £ + and £~ samples (rounded 
to the nearest integer), and expected statistical uncertainty on the absolute value of the measured 
lepton asymmetry, Ap B , for production-radiation tij events in pp collisions at Tevatron Run II, 
y/s = 2.0 TeV, for top quark mass m% = 178 GeV, summed over the two detectors with —0.1 < 
< 1.1 for the extra jet. K$ LO = for the tij rate. The upper block contains results for the 
dilepton sample, the middle block for lepton+jets sample, and the lower block for the combined 
samples. "Tight" cuts refers to those of Eq. |31 while "loose" refers to the possible increased 
acceptance as described in Sec. Ill Al The first four lines in each block represent improvements 
in detector and machine performance, while the last line represents uncertainty in the tij rate as 
discussed in Sec. lllll 
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We have also included a column for the value of the asymmetry, A £ p B , which is calculated 
directly from the cross sections and not the rounded numbers of leptons in the Table, and 
a column for the assumed K factor. For our very conservative baseline scenario with 
unlo/vlo — K = 1, the measurement is effectively consistent with zero. The very likely 
scenario of increased lepton acceptance and expected 5-tagging improvement, line 3, allows 
for a more interesting measurement. However, either increased luminosity or a fortuitous 
larger cross section, illustrated here with cr NLO /a LO = 2 corresponding to a smaller choice 
of factorization/renormalization scale, would be required to allow the statistical significance 
of the measurement to reach the level of 3 — 4 a. 

IV. ASYMMETRY IN TOP QUARK PAIR + JET EXCLUSIVE EVENTS 

We have seen that measuring the asymmetry in the inclusive ti sample may be doable but 
challenging, and that in the tij inclusive production-radiation sample it will be extremely 
difficult . The inclusive ti sample suffers from a smaller asymmetry, while for the inclusive tij 
sample statistics will be a larger problem. As seen in Sec. |H1 events without an additional 
hard jet should have a similar asymmetry to tij events, although of opposite sign. Yet this 
sample would be approximately 6 times the size of the inclusive tij sample, before any cuts 
are put on the jet rapidity. An obvious strategy is to experimentally attempt a jet veto on 
the ti inclusive sample to obtain an exclusive tiOj sample with better statistics than the 
inclusive tij sample. 

While definitely worth pursuing, this will not be straightforward. Recall that approx- 
imately as many events in inclusive ti production will have an additional hard jet from 
radiative top quark decay as from production radiation. One does not want to veto events 
with these jets. Certainly some fraction of the time D0 and CDF will be able to tell that the 
extra jet likely comes from radiative decay, based on the invariant mass of two or three jets 
reconstructing to a,W boson, or for a b jet and the extra jet plus either two jets or a lepton 
to have an invariant mass (or transverse mass) equal to m t . Determining how efficiently this 
can be done is beyond the scope of this paper, so we will ignore the complication of radiative 
top quark decays. Our goal in making estimates of how well one might measure A e FB in the 
exclusive sample is simply to highlight the beneficial features of the larger asymmetry and 
the greater statistics. 

Our procedure is to take the numbers of events expected in the inclusive ti sample based 
on the fraction of events that pass the cuts, as in Sec. |H1 and subtract from those forward 
and backward event numbers those from the tij sample. The obvious deficiency in this is 
that the inclusive asymmetry is known only from a NLO calculation (necessarily) which does 
not include decays, so the asymmetry shift from kinematical cuts effects cannot yet be taken 
into account. Thus, one should not put too much faith in the exact values predicted, as we 
expect future calculations to shift them noticeably. Our calculations, however, do represent 
the state of the art at this time. 

We show our numerical results in Table IIVI where the number of contributing events is 
approximately 18 times larger than in Table ITTT1 due to the larger cross section (x6) and 
the absence of a jet rapidity cut (x3). As before, we consider tight and loose levels of 
kinematic cuts, low and high Tevatron integrated luminosity, pessimistic and optimistic b- 
tagging efficiencies, and a possible 1-6-tag strategy for the dilepton sample only. Note that if 
the tij normalization is at the upper end of our cross section estimate, a factor of two larger 
than the conservative scale choice for the calculation, then the statistical significance of the 
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Sample 


/ Cdt [fbr 1 ] 


cuts 


£2b 


A FB 


Nf + N B 


N B +N e F 


..Hi 

K J 
NLO 


FB 
A FR 


signif. 


dilepton 


4 


tight 


0.20 


0.044 


176 


161 


1 


1.4 


0.7 a 


dilepton 


4 


loose 


0.20 


0.045 


339 


310 


1 


1.0 


1.0 cr 


dilepton 


4 


loose 


0.38 


0.045 


644 


589 


1 


0.75 


1.3 cr 


dilepton 


8 


loose 


0.38 


0.045 


1288 


1178 


1 


0.53 


1.9 a 


dilepton 


8 


loose 


0.62* 


0.048 


2891 


2644 


1 


0.35 


2.8 cr 


dilepton 


8 


loose 


0.38 


0.052 


1207 


1087 


2 


0.53 


1.9 a 


lepton+jets 


4 


tight 


0.20 


0.050 


502 


454 


1 


0.85 


1.2 a 


lepton+jets 


4 


loose 


0.20 


0.052 


764 


689 


1 


0.69 


1.4 a 


lepton+jets 


4 


loose 


0.38 


0.052 


1452 


1308 


1 


0.50 


2.0 cr 


lepton+jets 


8 


loose 


0.38 


0.052 


2903 


2616 


1 


0.35 


2.8 cr 


lepton+jets 


8 


loose 


0.38 


0.071 


2521 


2188 


2 


0.38 


2.6 a 


combined 


4 


tight 


0.20 


0.049 


678 


615 


1 


0.73 


1.4 Cr 


combined 


4 


loose 


0.20 


0.050 


1103 


999 


1 


0.57 


1.7 Cr 


combined 


4 


loose 


0.38 


0.050 


2096 


1897 


1 


0.42 


2.4 o- 


combined 


8 


loose 


0.38 


0.050 


4191 


3794 


1 


0.29 


3.4 o- 


combined 


8 


loose 


0.38* 


0.048 


5794 


5260 


1 


0.25 


4.0 Cr 


combined 


8 


loose 


0.38 


0.065 


3728 


3275 


2 


0.32 


3.2 cr 



TABLE IV: Numbers of forward and backward lepton events combining l + and £~ samples (rounded 
to the nearest integer), and expected statistical uncertainty on the absolute value of the measured 
lepton asymmetry, Ap B , for production-radiation- vetoed ttOj exclusive events in pp collisions at 
Tevatron Run II, = 2.0 TeV, for top quark mass rat = 178 GeV, summed over the two detectors. 
K% j LO = Z^iLo for the tfj rate The 

upper block contains results for the dilepton sample, the middle 
block for the lepton+jets sample, and the lower block for the combined samples. "Tight" cuts refers 
to those of Eq.|31 while "loose" refers to the possible increased acceptance as described in Sec. Ill Al 
The first four lines in each block represent improvements in detector and machine performance, 
while the last two line represents uncertainty in the ttj rate as discussed in Sec. II 111 The fifth line 
in the first and last blocks (the * entries) represent using a single 6-tag strategy for the dilepton 
sample as discussed in Sec. Ill Bl 

exclusive tiOj sample goes down slightly, due to the decreased sample size, which is partially 
compensated by the larger remaining asymmetry. As expected, the asymmetry in this 
sample is positive and somewhat larger in magnitude than the fully inclusive expectation of 
3.8%. While the statistics are better than indicated in Table ILTT1 these results still suggest 
that with the standard luminosity, the conservative cross section and loose cuts, only a 2.4 
a measurement is possible. An integrated luminosity greater than 4 fb -1 will be required 
for a more accurate determination of the asymmetry. 

V. DIFFERENTIAL ASYMMETRY IN TOP QUARK PAIR + 2 JET EVENTS 

The LO cross section for ttjj production at Tevatron Run II varies from 160 to 300 fb 
over the range of scale choices discussed earlier. We again calculate this with exact matrix 
elements from madgraph, choosing pr{j) > 20 GeV and < 2.0 as the limits in phase 

space integration for the extra jets. We further require AiL-j > 0.5 for the extra jet pair 
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FIG. 7: Normalized differential tijj cross section with respect to the higher-p^ extra jet pseudo- 
rapidity (dashed), overlaid with the total forward-backward positive lepton asymmetry (solid), in 
lepton+jets events with loose cuts as described in the text. 



to be separated from each other, to avoid the collinear singularity from gluon splitting. 
Compared to the 1-2 pb cross section for tij production, this is approximately a factor 5 
lower in statistics. Because of this, here we consider only the lepton+jets channel with the 
looser cuts and the improved b-tagging scenario. 

Similar to the plot in Sec. IIII Bl we show the normalized differential cross section with 
respect to the rapidity of the higher-pr of the two extra jets, as well as the differential 
asymmetry, in Fig. [TJ The asymmetry shows similar structure to that of Fig. El although here 
we use larger binning as it takes considerably more computing time to obtain the same level of 
statistics for smaller bins. The optimal region to use to obtain the best statistical significance 
for an asymmetry measurement is about —0.5 < \f](j)\ < 2.0, which includes approximately 
2/3 of the total rate. For 8 fb _1 per experiment and two combined experiments, with 
conservative cross section normalization and a double b-tagging efficiency of 38% Run II 
(summing the 2 detectors) could expect about 76 events with average asymmetry A l p B = 
—0.11 over this region. This would yield only about a 1.0 a measurement. If the cross section 
is really a factor of two larger (K = 2), this would still only be about a lAa measurement: 
interesting and worth pursuing, but difficult to achieve a useful level of precision. 



VI. COMPARISON WITH PARTON SHOWER MONTE CARLO 

Since Parton-shower Monte Carlo simulation packages, such as pythia 0], herwig 
and sherpa ^6| , use only LO matrix elements, which do not exhibit the interference struc- 
ture of the NLO matrix elements, they are not expected to reproduce the inclusive asym- 
metry seen at NLO. There is, however, the possibility that certain choices of the parameters 
that control the simulation of "color-coherence" in the showering/hadronization components 
of the Monte Carlo can lead to correlations that mimic the tij and tiOj asymmetries ob- 
served in the NLO perturbative results. Recall that, while the individual partons describing 
the short distance scattering carry color charges, the long distance initial- and final-state 
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hadrons are all color neutral. The requirement of reassembling at long distances into color- 
neutral states requires correlations between the partons. Complete fixed-order perturbative 
calculations account explicitly for the color- coherence, at fixed order, but Monte Carlo sim- 
ulations include color-flow correlations only approximately, based on various models. To 
illustrate this, we have studied ti dilepton production using pythia jl4j . 

We imposed the cuts of Eq.|3]on the b partons, leptons, and missing transverse momentum 
in Pythia-generated ti dilepton events. To identify extra radiation, we used the PYCELL 
subroutine with a cone size of R = 0.4. Any reconstructed jets were required to be AR > 0.4 
away from the b partons to be considered "extra radiation" . As the matrix elements used 
by pythia are LO, the lepton and antilepton distributions are the same, and the fully 
inclusive asymmetry is zero. The situation is more interesting when extra radiation is 
demanded or vetoed. In particular, correlations can be introduced by the angular structure 
of the all-orders QCD showering simulated in the Monte Carlo. For example, in pythia 
the parameters MSTJ(50) and MSTP(67) control the structure of the final state (MSTJ(50)) 
and initial state (MSTP(67)) showers. With the default values (the "on" value) for these 
two parameters, pythia adds initial- and final-state radiation to events based on color- 
flow information, while for the "off" value the showers develop independently of the color 
information and largely independently of the rest of the event. The color-flow constraint for 
the "on" case has the effect of restricting radiation to appear preferentially in the angular 
region between the p direction and the produced top quark direction, or between the p 
direction and the anti-top quark direction. Thus the presence of extra radiation tends to 
be correlated with larger angle scattering of the top quark, i.e., the top quark tends to 
be in the backward hemisphere when radiation is present. This correlation is intended to 
approximate the structure expected as a result of the "color coherence" in a full, all orders 
matrix element analysis with color singlet asymptotic states. Not only can this feature 
of pythia lead to a non-zero lepton forward-backward asymmetry in the presence of extra 
radiation (even when integrated over the rapidity of the radiation), but this asymmetry has 
the same sign and general correlation with the radiation as in the earlier tij perturbative 
analysis, as we now explore. Note that, by implication, the complementary event sample 
(with a veto on radiation) has an asymmetry of the opposite sign if the entire sample is to 
have a zero inclusive asymmetry. This means that the exclusive sample in pythia will have 
an asymmetry of the same sign as the tiOj exclusive sample. 

As in our earlier discussion of the matrix element calculations, CP invariance re- 
quires that the l~ distribution is always mirror-antisymmetric to that of the £ + about 
v{j) = 0; Ap> B (r) (j)) = —A FB (—i](j)). If the color-flow correlations in pythia 
are turned "off" (MSTP(67) and MSTJ(50) "off"), the final state has no knowledge of 
the C (non-invariant) structure of the initial state and (as in the earlier discussion) 
A l FB (r] (J) , off) = A FB (ji (J) , off) leading to a mirror-antisymmetric function of the jet ra- 
pidity A FB (rj (J) , off) = —A FB (—77 (J) , off). Thus the asymmetry need not vanish locally 
in 7] (j), but must yield zero when integrated over a symmetric interval of the jet rapidity 
( — Vo > > Vo)- The px correlation of Eq. [7| is such an integrated quantity and is 

exhibited in Fig. |HJ This figure provides a comparison of A FB between the pythia results 
with the color-flow correlation "off" (MSTP(67) and MSTJ(50) "off", the solid curve), PYTHIA 
results with the color-flow correlation "on" (the dashed curve) and the results of the NLO 
matrix element calculation (the dot-dashed curve), which we already considered in Fig. |21 
As expected, the "off" curve is consistent with zero, while the "on" pythia result and the 
NLO matrix element calculation exhibit comparable, negative values for A l FB (A FB has 
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FIG. 8: Asymmetry A^p B as a function of extra jet pt with pythia color-flow correlations turned 
off (solid) and on (dashed), as compared with the matrix element prediction (dot-dashed). 




FIG. 9: Asymmetry A l p B as a function of extra jet pseudorapidity rj(j) with pythia color-flow 
correlations turned "off' (solid) and "on" (dashed), as compared with the matrix element prediction 
(dot-dashed). 



equal magnitude but the opposite sign) with little j>t dependence. 

Similarly to our earlier discussion, the variation of the asymmetry as a function of the 
jet rapidity is a richer subject as displayed in Fig. As already argued the pythia 
"off" results show a mirror-antisymmetric function of the jet rapidity, A FB (j) (j) , off) = 
—A FB (—7] (j) , off), which integrates to zero over a symmetric rj (j) interval. In this case the 
"no radiation" exclusive sample will exhibit no asymmetry. Note that the "off" distribution 
in Fig. 01 is similar in shape to the gg subprocess results in Figs. El and but with the 
opposite sign. The sign of the asymmetry correlation can be understood by noting that 
in the color-flow correlation "off" pythia calculation the asymmetry arises predominantly 
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from cases where the extra jet comes from final state radiation (FSR) with the jet and 
the ti pair in the same hemisphere. This is to be contrasted with the gg subprocess 
NLO matrix element calculation where the extra jet is recoiling from the ti pair and tends 
to be in the opposite hemisphere. As indicated in Fig. El for the color-flow correlation 
"on" case (and the real matrix element result) the asymmetry knows about the C non- 
invariant structure of the initial state allowing a variety of features: A l FB (77 (J) , on) ^ 

A fb (v ti) , on), A FB (77 (j) » on ) ^ ~ A fb (~V ti) , on )> and 4b (V ti) , on ) y ields a non-zero 
net value when integrated over the jet rapidity. The fact that the net asymmetry for l + 
is negative (positive for t~) arises not from these symmetry considerations, but rather the 
details of the correlations. The corresponding "no radiation" complement has an asymmetry 
of the opposite, i.e., positive, sign (for £ + ). 

We can analyze the pythia calculations in somewhat more detail by separately turning 
"on" and "off" color-flow correlations, initial state radiation (ISR, as defined in pythia) and 
final state radiation (FSR, as defined in pythia). (Note that, due to interference effects 
there is no directly analogous analysis possible for the matrix element calculation.) The 
corresponding curves are exhibited in Fig. [TUlfor both £ + (solid) and £~ (dashed). (The 
case with all parameters set to "on" is the dashed cure in Fig. El and the case with all 
radiation "on" but color- flow correlations "off" is the solid curve in Fig. El) As suggested 
earlier, we see in Fig. a) and c) that with no color-flow correlations A FB (77 (j) , off) = 
A FB (77 (j) , off). On the other hand in Fig. El b) and d) with color-flow correlations present 
we have A l FB (77 (j) , on) 7^ A FB (77 (j) , on). From Fig. ITUla) we learn that ISR without color- 
flow correlations yields just a small (mirror-antisymmetric) asymmetry with the radiation 
in the opposite direction from the ti pair, while the addition of color-flow correlations in 
b) yields a much larger asymmetry with the t quark recoiling from (moving opposite to) 
the radiation (with the i approximately at rest) for forward radiation and the i moving 
forward (with the t approximately at rest) for backward radiation. It is this structure that 
is qualitatively similar to the qq subprocess matrix element result that dominates in Fig. |21 
Comparing Fig.lTUlc) and d) we see that the FSR contribution corresponds to the ti pair and 
the radiation all tending to move in the same direction, and that this mirror-antisymmetric 
structure is only slightly modified when color- flow correlations are present. Further, it is the 
FSR contribution that defines the qualitative shape of the combined result with no color-flow 
correlations present (the solid curve in Fig. EJ) . The combined result including color-flow 
correlations, the dashed curve in Fig. El has a shape that interpolates between Fig.HUJc) and 
Fig- Eld). However, since the FSR contribution in d) is close to mirror-antisymmetric, it 
is the ISR contribution in b) that dominates any symmetric integral over this distribution. 

The lesson here is that, even though pythia does not (cannot) exhibit the inclusive 
forward-backward asymmetry characteristic of NLO QCD, pythia does include parameters 
that can be set to produce events with much of the correlated asymmetry structure of 
NLO QCD. However, the agreement is only qualitative. In detail the asymmetries in the 
pythia event samples have neither the same magnitude nor the same dependence on the 
extra radiation as the matrix element prediction, i.e., the dashed and dot-dashed curves in 
Figs. |H1 and El are distinct. At least one underlying difference is that pythia includes FSR 
from the b quark in the top decay, which we do not include in the matrix element calculations. 
To the extent that the color-flow structure in pythia correctly reflects coherence effects in 
the higher order showering corrections, beyond NLO, the asymmetries exhibited by pythia 
are suggestive of what higher-ord er p erturbative analyses will yield. 

We have not studied herwig |l5( extensively except to confirm that the inclusive asym- 
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FIG. 10: Asymmetry A e FB as a function of extra jet pseudorapidity 77 (j) for £ + (solid) and t~ 
(dashed) for the following choices of pythia parameters: a) with color-flow correlations turned 
"off' , ISR "on" and FSR "off' ; b) with color-flow correlations turned "on" , ISR "on" and FSR 
"off'; c) with color-flow correlations turned "off', ISR "off' and FSR "on"; d) with color-flow 
correlations turned "on", ISR "off' and FSR "on". 

metry is not present. It may or may not have correlated asymmetries similar to those 
observed with pythia arising from color-flow related constraints on extra radiation from 
showers. We have performed no analogous analyses with sherpa ^|. However, future 
studies using these generators should be aware that these correlated asymmetries may be 
present. In any case, it is clear that the experimental study of asymmetries, both inclusive 
and exclusive, can yield information about the coherence structure of QCD showers. To 
the extent that top quark physics is a background in studies of physics beyond the Standard 
Model, we must come to understand quantitatively these correlations and asymmetries in 
Monte Carlo simulations. 

VII. CONCLUSIONS 

Top quark pair production at a pp collider exhibits a forward-backward asymmetry due to 
higher-order short-distance QCD effects. This has been known for almost two decades, but 
only at an abstract theoretical level. Contributions to the asymmetry arise both from virtual 
corrections to qq — > tt, due to interference between the color-singlet box (4-pt.) and Born 
terms, and from real gluon emission in the same subprocess, due to interference of emission 
between the initial and final states (although one cannot think of it physically so simply, 
because those and other diagrams for emission are gauge- related) . The two contributions 
are opposite in sign and of similar magnitude, leading to a partial cancellation which leaves 
a smaller asymmetry in the inclusive sample of the same sign as the virtual piece. 
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We have explored the asymmetry here with an eye toward experimental measurement. 
There are several major new points. First, experimentally the most reliable avenue for 
measuring the asymmetry is likely via the charged leptons, where both the direction and 
the sign of the charge can be determined with confidence. Unfortunately, we cannot test 
this approach in detail for an inclusive ti sample without a (currently unavailable) full NLO 
calculation including the t decays. For the tij sample (requiring extra radiation), where the 
asymmetry appears at lowest non-trivial order, our results do confirm that the asymmetry 
at the top quark level is reflected in the top quark decay products, although with different 
magnitude. This change in the magnitude is due to the top quarks' spins leading to decay 
products being emitted preferentially in certain directions, which can move them from one 
hemisphere to another. Secondly, we have pointed out that the asymmetry of the leptons 
in the tij sample is correlated with the angular distribution of additional jets in the event, 
but not their transverse momentum. Unfortunately, the real-emission asymmetries, which 
we term the tij and tijj components, are presently known only at LO. These last aspects of 
the short- distance asymmetry were not previously noticed in the literature. Finally we have 
studied the likely numbers of events at the Tevatron for the fully inclusive ti sample, where we 
assume the leptonic asymmetry is equal to the top-quark asymmetry, the radiation-required 
ttj sample, and the jet- veto-based exclusive tiOj sample. In all cases statistics will be an 
issue at the Tevatron, requiring a multiple channel analysis with (hopefully) outstanding 
accelerator and detector performance. 

We have also shown that some models for color-flow correlations in QCD showers can yield 
correlated asymmetries as well, as exemplified in the parton shower Monte Carlo pythia[i4^. 
While the lepton asymmetry always vanishes in the totally inclusive PYTHlA-generated sam- 
ple, when the relevant shower parameters are set (to the default values) to exhibit a specific 
model of color-flow structure in the parton showers, a negative top quark asymmetry arises in 
the tij component with a positive asymmetry in the corresponding exclusive Oj component. 
Both the overall magnitude of the asymmetry and the differential asymmetry with respect 
to the additional hard jet angular distribution are qualitatively similar to but different in 
detail from those predicted by the perturbative NLO results. 

We again emphasize that we have tried to explore only the general structure of the ti 
asymmetry signal and that, due to a number of limitations in the current theoretical tools, 
our results should not be regarded as precise predictions. We review here the major sources 
of those uncertainties/questions, which should be considered as areas that theorists need to 
address in the near future: 

• How does the inclusive asymmetry change once kinematic cuts are imposed on the 
decay products? No NLO calculation exists to address this question. 

■ What is the overall tij normalization, which is known presently only at LO? Ad- 
ditionally, we need to evaluate the NLO rate including decays so that cuts may be 
imposed. 

■ What is the effect of radiative top quark decays on experimental efficiencies for cor- 
rectly selecting the extra jet in the inclusive production-radiation tij sample? We need 
this to be able to isolate the maximally-asymmetric region in r)(j). 

■ We need to properly merge the tij matrix elements into the parton-shower Monte 
Carlo environment in order to be able to make complete predictions. 
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While the experimental measurement of the asymmetries described here will likely be a 
challenge at the Tevatron due to limited statistics, it is still an extremely important goal. 
Since top quarks may well have a unique connection to new physics, and since they play a 
large role as a background in many new physics searches, Tevatron experimentalists should 
strive to determine top quark properties as accurately as possible. Further, the analysis 
of the correlation between the asymmetry and extra radiation may offer a nearly unique 
opportunity to elucidate detailed properties of long-distance, all-orders QCD, i.e., the color- 
flow structure of the parton shower. While the measurement of A £ FB is unlikely to be a direct 
probe for new physics, unless the new physics contributions to the asymmetry are extremely 
large or of opposite sign to the SM, this measurement still provides a good example of the 
subtle behavior of SM particles in unexpected places. We should understand the SM as well 
as possible, from both theory and data, before moving on to see what lies beyond. 
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